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This  report  was  or  1 g I na I  I y  prepared  and  distributed 

sr,ss Ms;.';;:  s: 

for  a  formal  reference. 


This  UC1D  contains  the  original  text  and  re“ ' ^ed . “u 1 E JJe 
f..  11.  >est  cases  S...r.j  c  anges  were  «d.  ^ 

?o  ?o«d°o:?  'for1  00  cm  P?STs  »ersus  X  are  no.  on  a 

to  facilitate  exlrac  Ion  of  sma  I  X 
results.  Snapshot  output  Is  now  presented  at  1.0,  W.U, 
20.0,  50.0.  and  100.0  cm.  In  X. 


The  results  for  lest  case  B  as  original ly  distributed 
were  n  e^ror  as  pointed  out  by  Dr.  Keith  BruecUne-  The  X 
step  employed  .as  loo  large  to  pcoperlv  describe  the 
avalanche  Ionization.  Hence.  In  this  report 
been  run  with  reduced  X  steps. 


al  I  cases  have 
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'  V  At  the  high  current  beam  propagation  workshop  (Ref.  1 )  It  was 
realized  that  fleldsojjer  comparisons  betwssn  working  groups  should  be 
made  for  booms  with  oxpandod  heads  such  as  develop  by  blowoff  and 
erosion  during  propagation.  To  facilitate  comparisons  of  the 
fieldsolvers  while  removing  the  differences  due  to  evolution  of  the 

beam  during  propagation^  standard  test  cases  with  imposed  beam  head 

0»‘-  '  < 

profiles  twnr» - been^  suggested.  Hence,  unambiguous  field  and 

conductivity  generation  calculation  comparisons  can  be  made. 


The  field  equations  solved  for  this  report  are  derived  elsewhere* 
(Ref.  2).  The  plasma  generation  model  includes  a  direct  deposition 
term,  an  avalanche  ionization  term  based  on  experimental  data  (Ref.-? 
3),  and  a  recomb i na t i on  term.  Plasma  density  is  reloted  to  a  scalar 
conductivity  through  the  momentum  transfer  collision  f  requency  . -  These ... 
models  have  been  used  by  the  Beam  Research  Group  at  LLNL  in  the 
EMPULSE  I  (Ref.  4).  EMPULSE  II  (Ref.  5),  RINGBEARER  I  (Ref.  6),  and 
RINGBEARER  II  (Ref.  7)  codes.  O-The  actual  coded  fieldsolver  used  for 
the  results  in  this  report  was  extracted  from  the  EMPULSE  II  code. 
The  (cod^i  PHOENIX^  was  constructed  to  run  this  fieldsolver  and  generate 
the  appropriate  diagnostic  output. 


Standard  test  cases  are  proposed  where  the  beam  current  density 
is  prescribed  at  all  points  in  space.  This  is  accomplished  using  a 
simple  parameterized  profile  which  can  represent  a  cylindrical  beam  or 
one  with  an  eroded  nose.  Four  test  cases  are  proposed;  PHOENIX 


results  for  these  cases  are 
quantities  are  in  cgs  units; 


presented .  ■*—  In 

>eom  ” temporary 


this  report  almost  al I 
parameters  such  as  rise 


time  (length)  are  always  given  in  centimeters.  \ 
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Ms^lLLLds 

Coordinate  System  -  Z.  X.  r.  8: 


Z  -  distance  of  propagation  of  a  beam  segment  -  since  beams  are 
not  propagated  in  the  test  cases  Z  =  0. 

X  -  distance  back  from  the  beam  head. 

r.0  -  transverse  coordinates  in  the  cylindrical  geometry  - 

since  only  monopole  fields  are  calculated  all  quantities  are 
0  independent. 

Field  Equations: 


13  3  An  dd  (1) 

- r  —  &+<p)  = - — ) 

r  3r  3r  c  3X 


1  3  dl4  1  3  A no  dip  (2) 

r  3r  3r3X  r  dr  c  dr 


Boundary  Conditions: 


3< d  dip 

— ( r=0)  =  — (r=0)  =  0 
3r  3r 


(3) 


jd(  r=b)  =  <p(  r=b)  =  0 


(4) 
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j4  and  <p  are  the  usual  monopole  field  potentials,  a  Is  the  plasma 
conductivity,  and  J8Z  is  the  Z  directed  component  of  the  beam  current. 
The  parameter  b  is  the  radius  at  which  there  Is  assumed  to  be  a  metal 
wall.  A  term  of  the  form  (d*j4/dX*)  has  been  dropped  In  equation  (1). 
except  in  the  RINGBEARER  I  code  where  this  term  is  retained.  The 
RINGBEARER  I  code  uses  a  field  equation  obtained  by  combining 
equations  (1)  and  (2)  and  using  r2  as  the  independent  variable  in  the 
radial  direction.  In  this  case  the  boundory  condition  on  jA  at  r  =  0 
is  modified  (Ref.  6). 

Conductivity  Calculation: 

The  simplified  air  chemistry  models  and  optional  features 
originated  by  E.  Lea  (Ref.  4,5)  are  described  in  the  following 
section.  Several  minor  changes  have  been  made  In  the  recombination 
and  collision  frequency  models  to  allow  the  extraction  of  a  single 
electron  temperature  parameter  and  to  simplify  the  analytic  forms. 
These  changes  will  not  affect  computed  conductivities.  The  field 
dependent  recombination  and  collision  frequency  models  are  included 
for  completeness  since  they  have  been  used  in  the  past.  However, 
current  understanding  of  the  air  chemistry  indicates  better  models  are 
needed  which  retain  the  simplicity  and  utility  of  these  models  but 
with  increased  accuracy.  Hence,  we  are  not  advocating  the  use  of 
these  models.  In  the  standard  cases  reported  herein  recombination  is 
omitted  and  the  high  electron  temperature  limit  of  the  collision 
frequency  is  used. 

The  generation  of  electron  density,  n#,  Is  given  by: 
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Hare  a  =  4.803  10-10  esu  and  c  =  2.9979  10,fl  cm/s.  Typically, 


d£/dX  =2.52  MaV  cm2/gm  (7) 


...  — — — — m  poll  II I|JI •mPOTmifppiHMMWlI 


E2 

S  =  =- 


(12) 


where  E  is  the  magnitude  of  the  electric  field  in  statvolts/cm  and  p 
is  the  gas  density  normalized  to  that  at  273°K  and  one  atmosphere. 

For  the  computation  of  the  recombination  rate,  a,  and  the 
momentum  transfer  collision  frequency,  um,  an  electron  temperature, 
T(,  is  calculated  from  E/p.  Temperature  is  assumed  to  vary  linearly 
with  E2/p2  due  to  ohmic  dissapation;  when  E  =  0  the  electron 
temperature  is  the  gas  temperature. 


T.(eV)  =  T  (eV)  +  *  (E2/p2) 


(13) 


The  value  of  k  is  determined  in  slightly  different  manners  for  the 
recombination  and  collision  frequency  calculations.  The  intent  in 
computing  T#  is  as  a  parameter  for  computing  i/m(E/p)  and  a(E/p);  and 
T#  is  not  intended  to  accurately  indicate  the  electron  temperature  for 
all  E/p.  Since  the  formula  for  asymptotes  above  T,  «<  2  eV  and 
above  this  temperature  a  is  quite  small  equation  (13)  is  an  adequate 
approx imat ion. 


Recombination,  when  included.  Is  a  function  of  E/p: 

2.10  10-7  2.10  10“7 


a(cm3/s)  = 

(T,(eV)/.02586)-7  (T#(°K)/300. ) • 7 

*  =  4.6774  10-3 


(14) 


This  value  of  k  was  determined  by  requiring  that  T#  =  2/3  eV  when  = 
.02586  and  £2/p2  =  137.0.  This  recomb  I  not  Ion  rote  Is  due  to  the 
presence  of  02*  and  was  suggested  by  Bob  Johnston  (Ref.  8).  For  the 
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standard  casa  results  thus  far  computed  recombinat Ion  has  been 
omitted. 


Assuming  a  scalar  conductivity,  a  is  related  to  n,  by: 

n.»2  (15) 

m  =  9.11  10-2®  grams  is  the  electron  moss;  vm  is  the  momentum  transfer 
collision  frequency: 


The  momentum  transfer  collision  frequency  is  modelled  as  a 
funct i on  of  E/p  as : 


_  A  +  T 

v m(s"1 )  =  A. 86  1012  p  - 

8  +  T# 

*  =  4.8662  10-3 
A  =  2.35835  10~3 
B  =  .561026 


(16) 


This  value  of  k  was  determined  by  requiring  T#  =  2/3  eV  when  Tg  =  0 
and  E 2/p2  =  137.0.  This  formula  was  determined  by  requiring  agreement 
with  tables  (Ref.  9)  at  T#  =  .0235  eV,  2/3  eV,  and  as  T#  ->  °°.  In  the 
standard  case  results  the  high  temperature  limit  (T#  >  2  eV)  was  used; 
name  I y , 

i'm(s"')  =  ‘♦•86  10<2  p  (17) 


When  recombination  is  ignored  and  vm  is  constant  n#  —  a  and  the 
conductivity  con  be  calculated  directly  using: 


-  8  - 


Beam  Prof i I e : 


I8(*)  =  I80  »anh(X/LR) 

IB0  -  Beam  current  at  X  =  « 

Lr  -  rise  length  in  centimeters 

JB(r.X)  =  J„(r/R(X))  IB(X)/(rrR2(X) ) 

R(X)  =  ROMIN  +  { R0MAX-R0MIN)/2 . 0  ( t-tanh((X-XLOC)/XWDTH) 

ROMIN  -  Beam  radius  as  X->°° 

ROMAX  -  Beam  radius  as  X->-°° 

XLOC  -  Location  ot  "midpoint"  of  neck 
XWDTH  -  Width  of  neck  region 

Gaussian  -  J(r/R)  = 

Bennett  -  J(r/R)  =  1/(l+r2/R-)2 

Gas:  760  torr  (1  atm)  air  and  76  to rr  (.1  atm)  air 

i/  constant 

m 

a  =  0 
Numer ica  (  : 

XMAX  =  200.0  cm. 

Output: 

Snapshot  profiles  in  r  of: 

JB’  JNET'  ~JRET:  u t •  c/" ' • 

:  V  •  ^  r  >  ^  z 

at:  X  =  10. .  20. ,  40. ,  100. ,  200.  cm. 


History  profiles  in  X  of: 

IB;  Rx :  IB .  Inet*  ^ret*  ^efe*  ^z^r=®)*  <7(r=®)>  */j(r=0) 

Contour  plots  in  X,r  of; 
jtf(X.r);  i/  ( ( X  ,  r ) 
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IV.  Standard  Cases 

10  kA,  Standard  Case 

A  10  kA  Bennett  beam  into  full  atmosphere  at  injection  - 
beam  is  cylindrical  with  no  eroded  head. 

IB0  =  10000.0  A  R0MIN  =  .5  cm 

'-RISE  =  '0-0  cm  ROMAX  =  .5  cm 

y  =  100.0  XLOC  =  15.0  cm 

p  =  760  torr  XWDTH  =  15.0  cm 

10  kA,  Reduced  Pressure 

A  10  kA  Bennett  beam  into  one  tenth  atmosphere;  at  injection  - 
beam  is  cylindrical  with  no  eroded  head. 

IB0  =  10000.0  A  R0MIN  =  .5  cm 

Lrise  =  10.0  cm  ROMAX  =  .5  cm 

y  =100.0  XLOC  =  15.0  cm 

p  =76  torr  XWDTH  =  15.0  cm 

100  kA,  Broad  Nose 

A  100  kA  Bennett  beam  with  an  imposed  eroded  head  but  with  a 
broad  nose. 

IB0  =  100000.0  A  R0MIN  =  . 2  cm 
'-RISE  =  cm  ROMAX  =  40.0  cm 

y  =100.0  XLOC  =  15.0  cm 

p  =  760  torr  XWDTH  =  15.0  cm 


100  kA,  Narrow  Nose 

A  100  kA  Bennett  beam  with  an  eroded  head  with  radius  parameters 
fitted  to  an  SAI  run. 

IB0  =  100000.0  A  R0MIN  =  .5  cm 

LRISE  =  300  cm  ROMAX  =  100.0  cm 

y  =100.0  XLOC  =  20.0  cm 

p  =  760  torr  XWDTH  =  1 .8  cm 
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V.  PHOENIX  Output 


F  rome 

1.  Input 

parameters  are  defined  as 

foil ows : 

C- 

C- 

r_ 

NAMCCISr 

INPUT  FOR  THE 

PHOENIX  FIELDSOLVER 

c- 

NAME 

UNITS 

DEFAULT 

COMMENTS 

L.- 

c- 

BEAM  PHYSICAL  PARAMETERS 

c- 

IBO 

AMPS 

10000. .0 

BEAM  CURRENT 

c- 

LR 

CM 

30.0 

BEAM  RISE  LENGTH 

c- 

GAMMAO 

100.0 

BEAM  ENERGY  PARAMETER 

c- 

IPROFILE 

1 

=1  FOR  GAUSSIAN  PROFILE 

c- 

=2  FOR  BENNETT  PROFILE 

c- 

ROMIN 

CM 

.5 

MINIMUM  RADIUS  AT  LARGE  X 

c- 

ROMAX 

CM 

40.0 

MAXIMUM  RADIUS  AT  LARGE  -X 

c- 

XLOC 

CM 

15.0 

LOCATION  OF  NECK  REGION 

c- 

XWDTH 

CM 

15.0 

WIDTH  OF  NECK  REGION 

c- 

GAS  PARAMETERS 

C- 

ENOO 

CM»*-3 

1  .  E4 

INITIAL  ELECTRON  DENSITY 

c- 

KSCATTER 

1 . 699 1 57E3 

SCATTERING  PARAMETER 

c- 

4(E»*2/MC»*2)(CN/MC*»3) 

c- 

KBEAM 

6.69 

DIRECT  DEPOSITION  PARAMETER 

c- 

( DE/DX  )  •  ROE/(  DE  L  T  AE  •  E  •  C  ) 

c- 

LAMBDARO 

CM 

3 . 1 3336E4 

RADIATION  LENGTH 

c- 

IFROHAT 

1 

-1  FOR  NO  CHANNEL 

c- 

r _ 

RHOHATO 

1 .00 

OENSITY.  VARY  TO  VARY  PRESSURE 

t- 

c- 

R  GRID 

c- 

c- 

NR 

NRMX 

NUMBER  OF  R  GRID  POINTS 

c- 

B 

CM 

10.0 

MAXIMUM  RADIUS  FOR  FIELDS 

c- 

NRLNSTEP 

40 

=0  PURELY  EXPONENTIAL  R  GRID 

c- 

=NR-1  LINEAR  GRID 

c- 

> 1 ,  <  NR-1 i  HYBRID  GRID 

c- 

RO 

1  .0 

ORIGINAL  BEAM  RADIUS 

c- 

RLNFAC 

1  .0 

FRACTION  OF  RO  TO  BE  LINEAR 

c- 

X  GRID 

c- 

c- 

NX 

NXMX 

NUMBER  OF  X  GRID  POINTS 

c- 

NDX 

1 

NUMBER  OF  DIFFERENT  SIZE  STEPS 

c~ 

DXVALUES 

CM 

.2 

ARRAY  OF  X  STEP  SIZES 

c- 

X LIMITS 

ARRAY  OF  SWITCHING  PLACES 

C- 
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Frames  2-11.  Snapshots  versus  r  at  several  X. 


Current  densities  are  J8,  JRrT=cE2,  and  JNet=jb+jRET •  Since 
usually  JRET<0,  -JRET  is  plotted. 

c/i/.  is  plotted  to  indicate  the  efolding  distance  in 
centimeters  for  the  computed  avalanche  ionization  rate. 

Frames  12.13  Prescribed  beam  current  and  radius  profiles. 

Frame  14  Neutral izat ion  fractions,  the  plotted  quantities  are: 

"B"  ->  IB(X)/It0 

"R”  ->  -Irej/Ibo  =  /ob  2?rrdr/1 B0 

"N"  ->  (Ib+iret)/ibo  =  inet/ibo 


'll"  ->  IEff/ibo  =  /ob 


er  drf  2nrdr  JR 
ymc2  dr  IB 


IEFr  is  the  effective  current  providing  the  beam  pinch  force; 
if.  in  the  magnetic  regime,  IREt  =  ®  (and  JRET=0)  then  IEFf.  = 
IB:  if  JRE r  has  the  identical  profile  to  JB  then  IEFF  ~  INEf 


Frame s  15,16.17.  X  histories  of  the  on-axis  (r=0)  values  of  Ez,  o, 
and  i/,  . 


Frame  18.  Contour  plot  of ^rf(X.r),  contours  are  at  . 1-.9 ^maK • 

Frame  19.  Contour  plot  of  i/,(X,r),  contours  logarithmically  spaced. 
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NUIOUT  RECOMBINATION  RATE  ALPHA 


TIME,  DATE: 

13:52:33 
U  03/1 1/81 


< 

x 

a. 

-j 

< 


Figure  4.  Plot  of  recombination  rate,  a,  versus  E/p.  E 
field  values  given  on  the  horizontal  axis  are  for 

p  =  1.0. 
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ENDTEXT 


PHOENIX 


BEAM  Pn»  AMETERS 

IBO(AMPS) 

=  1 . 000E+04 

LR(CM) 

1 . 000E+01 

GAMMAO 

-  1 . 000E+02 

IPROF 1 LE 

2 

ROMIN(CM) 

=  S.OOOE-OI 

romax(cm) 

5 . OOOE-OI 

XLOC(CM) 

-  1  .  500E  +0 1 

XWDTH(CM) 

1 . 500E+01 

1 S 1 GMOD 

=  0 

GAS  PARAMETERS 

ENOO(CM*  * 

3)  -  1  .OOOEf-04 

KSCATTER 

1  . 699E+03 

KBEAM 

:  6.690E+00 

LAMBDARO 

- 

3. 133E+04 

IERHOHAT 

) 

RHOHATO 

r: 

1 .OOOE+OO 

numerical 

parameters 

NR 

=  220 

B(  CM) 

- 

1 . 000E+01 

nrlnstep 

;  40 

RO(CM) 

1 .OOOE+OO 

Rl  Nf  AC 

-  1 . OOOE+OO 

NX 

- 

2500 

NDX 

-  1 

DXVALUES( t ) 

CM= 

4 . 000E-02 

CASE  :  A  FRAME :  I 
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J(ESU) 


Snapshot  K 


IO.OO,  *  STEP  WO 


»(CM) 


SNAPSHOT  * 


tO  00.  X  STEP  NO 


J(£SU) 


SNAPSHOT  *  =  2000,  X  STEP  NO 


t*T2  9  ©  <n  © 

Hot) 


Hot) 


R«S)  "  ~  »(CM) 


SNAPSHOT  v 


20  00.  X  STEP  NO 


PHOEN  I  X 


FIELDSOLVER  WITH  FIXED  JB(R.X) 
PLOT  VERSUS  X 


CASE  A  f«AMf 


*5 


X  (CM) 


PHOENIX 


TIME.  DATE 
13:36 : 58 
09/22/8 1 


l BEAM  = 

9 . 993E+03 
LRISE  = 

I . OOOE+O I 
ROM  IN  = 

5 . OOOE-OI 
ROMAX  = 
5.000E-01 
XlOC  = 

1  500E+01 
XWDTH  - 
1  500E+01 


CASE  A  F  RAM! 


PHOEN  I  X 


TIME,  DATE 
13:36:58 
09/22/8 f 


I  BEAM  - 

9  99Jf+0.5 
LR»SE  ~ 

t .oooe+or 

ROM  IN  - 
5.000E-01 
ROMAX  s 
5. OOOE-OI 
XLOC  = 

1  500E+-01 
XWDTH  = 

1  500E+01 


CASE  A  FRAME 


FIELDSOLVER  WITH  FIXED  JB(R,X) 
PLOT  VERSUS  X 


FIElOSOLVER  WITH  FIXED  JB( R  ,  X  ) 
PLOT  VERSUS  X 


X(CM) 


-  26  - 


i 


PHOEN  I  X 


FIELDSOLVER  WITH  FIXED  JB(R, 
PLOT  VERSUS  X 


TIME,  DAT l 
15:56:58 
09/22/81 


'BEAM  =  ' 

9  9951*03 
LRISE  s 
1  OOOE+OI 
ROM  IN  = 
5.000E-01 
HOUAX  - 
5 . 0001 -01 
XLOC  = 

1  . 500E  +01 

xwoth  = 

1  500E+01 


case  »  (»auc 


PHOENIX 


TIME,  DAI  E 
1 3 : 56  58 
09/22/81 


I  BE  AM 
9  995E+0J 
LRISE  = 

I  OOOE+OI 
ROM  IN  = 

5  OOOE-01 
ROMAX  = 

5  OOOE-OI 
XLOC  •- 
1  500E+01 
XWOTH  — 

1  500E+01 


X(CM) 

le 


FIELDSOLVER  WITH  FIXED  JB(R,X) 
PLOT  VERSUS  X 


CASE  A  ISAM!  17 


X(CM) 


£♦02 


F I E LDSOLVER  WITH  FIXED  J0(R,X) 
CONTOUR  PLOT  OF  AO(R,X) 


X(CM) 


TIElDSOLVER  WITH  FIXED  JB(R,X) 
CONTOUR  PLOT  OF  NUSUBl(R.X) 

I -r- ■■ -rr-n Tl 


PHOEN  I  X 


BEAM  PARAMETERS 

1  BO (AMPS) 

— 

1 . 000E+04 

LR(CM) 

= 

1 . 000E+01 

GAMMAO 

= 

1 . 000E+02 

1  PROF  1 LE 

= 

2 

ROMIN(CM) 

= 

5 . OOOE-OI 

ROMAX(CM) 

= 

5.000E-01 

XLOC(CM) 

- 

1 . 500E+01 

XWDTH(CM) 

= 

1 . 500E+01 

1 S 1 GMOD 

= 

0 

GAS  PARAMETERS 

ENOO(CM*  *- 

3)  - 

1  . 000E+04 

KSCATTER 

— 

1  . 699E+03 

KBEAM 

r 

6 . 690E+00 

LAMBDARO 

- 

3. 133E+04 

1  FRHOHAT 

1 

RHOHATO 

— 

J . OOOE-OI 

NUMERICAL 

PARAMETERS 

NR 

' 

220 

B(CM) 

— 

1  . 000E+01 

NRlNSTEP 

40 

RO(CM) 

= 

1 .OOOE+OO 

RLNF AC 

- 

1 . OOOE+OO 

NX 

= 

2500 

NDX 

- 

1 

DXVALUES( 1 ) 

CM- 

4 . 0001-02 

CASE .  B  FRAME :  1 


-  29  - 


/ 


PHOEN  I  X 


F I ELDSOLVER  WITH  FIXED  JB(R,X) 
PLOT  VERSUS  X 


TIME,  DATE 
13:52:50 
09/22/81 


I  BEAM  = 

9 993E+03 
LRISE  = 

1 . 000E+01 
ROM  IN  = 

5 . OOOE-Ol 
ROMAX  = 
5.000E-01 
XLOC  = 

1 . 500E+01 
XWDtH  = 

1  500E+01 


X(CM) 

CASE  6  f  RAMC  1 2 


PHOENIX 


TIME.  DAU 
13:52-50 
09/22/81 


I  BEAM  - 
9  993E  *03 
1  R  i  SE  - 
1  OOOE^OI 
ROMIN  = 

5  OOOE-Ol 
ROMAX  - 
5.000E-01 
XLOC  = 

1  . 500E+01 
XWDTH  - 
1  500E+01 


FIELDSOLVER  WITH  FIXED  JB(R,X) 
PLOT  VERSUS  X 


CASE  8  f Rami  13 


X  (CM) 


PHOENIX 


FIELDSOLVER  WITH  FIXED  JB(R,X) 
PLOT  VERSUS  X 


TIME,  DATE 
1  3 : 52 : SO 
09/22/81 


I  BEAM  = 

9 . 993E+03 
LR I SE  = 

1 .OOOE+OI 
ROM  IN  = 

S.OOOE-OI  £ 
ROMAX  =  £ 

5  000E-01  3 
XLOC  = 

I . 500E+0I 
XWDTH  — 

1 . SOOE+Ot 


CASE  8  fRAUE  14 


X  ( CM) 


PHOENIX  FIELDSOLVER  WITH  FIXED  J B ( R 

PLOT  VERSUS  X 


TIME,  DATE 
13:52:50 
09/22/81 


I  BEAM  - 
9  993E+03 
[RISE  = 

1  .0001*01 
ROM  IN  - 

5.000E-01  o 
ROMAX  =  ^ 

5.000E-01  ^ 
XLOC  ,xi 

t  50oe*oi 

XWOTH  - 
1 . 500E+01 


CASE  0  f  RAMI 


5 


X(CM) 


20+ 1  ■  rp - ■ - ■ - ■ - ■ - < - 1 - > - - - 1  TO*  3 


CASE  8  f  RAM(  1 8 


PHOEN  I  X 


TIME,  DATE 
13:52:50 
09/22/81 

<«UH  • 

i«i«  « 

■  am*«t 


FIELDSOLVER  WITH  FIXED  JB(R,X) 
CONTOUR  PLOT  OF  NUSUBl(R.X) 


PHOEN  I  X 


BEAM  PARAMETERS 


1 BO( AMPS) 

-  1  . 000E+05 

LR(CM) 

= 

3  .  OOOE+Ol 

GAMMAO 

=■  1  .0001+02 

1  PROF  1 LE 

= 

2 

ROM ) N( CM) 

-  2  .  OOOE-OI 

ROMAX(CM) 

- 

4 . OOOE+O 1 

XLOC(CM) 

•-  1.500E+01 

XWDTH(CM) 

= 

1 . 500E+01 

1 S 1 GMOD 

0 

GAS  PARAMETERS 

l  NOO( CM*  * 

-3)  1.000E+04 

kscatter 

1  699E+03 

KBf.  AM 

-  6.690E+00 

lambdaro 

= 

3.  133E+04 

1 ERHOHAT 

t 

RHOHATO 

= 

1 .OOOE+OO 

numerical 

PARAMETERS 

NR 

220 

0(CM) 

— 

1  OOOE+Ol 

NRLNSTFP 

40 

RO(CM) 

= 

1  OOOE+OO 

RLNI  AC 

1 . OOOE+OO 

NX 

= 

2500 

NDX 

1 

OXVALUCS( 1 ) 

CM— 

4 .0001-02 

CASE  C  FRAME .  < 


-  39  - 


T 


T 


PHOENIX 


E I ELDSOLVER  WITH  FIXED  JB(R, 
PLOT  VERSUS  X 


CAM  C  ’««*  <5 


X(CM) 


£♦02 


PHOENIX 


fieldsolver  with  fixed  jb(r, 

PLOT  VERSUS  X 


TIME ,  DATE 
13:57:15 
09/22/81 


IBEAM  = 
9.993E*04 
LRISE  = 

3 . OOOE  *01 
ROM  IN  = 
2.000 £  01 
ROMAX  r 
4  OOOE *01 
XLOC  = 

1 . 500E+01 
XWDTH  = 

1  500E  *01 


o 

n 

a 

ta 

i/> 


X(CM) 

CASE  C  ER**  16 


PHOENIX  FIELDSOLVER  WITH  FIXED  JB(R, 

PLOT  VERSUS  X 


TIME,  DAH 
13:57:15 
09/22/81 


IBEAM  - 
9  993E+04 
IRISE  = 

3  OOOE  *0 1 
ROM  IN  = 

2  OOOE -01 
ROMAX  = 

4  OOOE *01 
XLOC  = 

1  500E*01 
XWDTH  = 

1 . 500E  *01 


o 

n 

O: 


CASE  C  rt/MME 


X(CM) 


t+02 


PHOENIX 


E l ElDSOLVER  WITH  FIXED  JB(R,X) 
CONTOUR  PLOT  OF  AO(R,X) 


TIME.  DATE 
13:57:15 


PHOEN  I  X 


BEAM  PARAMETERS 

1  BO (AMPS) 

r: 

1 . 000E+05 

LR(CM) 

3 . 000E+01 

GAMMAO 

- 

1  .  OOOE+02 

1  PROF  HE 

= 

2 

ROMIN(CM) 

5 . OOOE-O 1 

ROMAX(CM) 

1 . 000E+02 

xloc(cm) 

2  . OOOE+OI 

XWDTH(CM) 

1 . 800E+01 

1 S 1 GMOD 

- 

0 

GAS  PARAMETERS 

ENOO(CM*  *- 

-3)  - 

1 .0OOE+O4 

KSCATTER 

1 . 699E+03 

KBEAM 

- 

6.690E+00 

LAMBDARO 

3. 133E+04 

1  ERHOHAT 

= 

1 

RHOHATO 

= 

1 .OOOE+OO 

NUMERICAL 

PARAMETERS 

NR 

- 

220 

B(CM) 

zz 

1 . 000E+01 

nrlnstep 

- 

40 

RO(CM) 

= 

1 .OOOE+OO 

RLNFAC 

1 . OOOE+OO 

NX 

= 

2500 

NDX 

= 

1 

DXVALUES( 1 ) 

o 

r 

n 

4 . 000E-02 

CASE :  D  FRAME :  1 


-  49  - 


PHOEN  I  X 


FIELDSOLVER  WITH  FIXED  JB(R.X) 
PLOT  VERSUS  X 


TIME.  OAfE 
14:04:23 
09/22/81 


I  SEAM  = 
9.993E+04 
1R1SE  = 

3  OOOE+OI 
ROM  IN  = 
5.000E-01 
ROMAX  = 

I  OOOE+02 
XLOC  = 
2.000E+01 

xwoth  = 

1  800E+01 


1/1 

LU 


X(CM) 


CXSC  0  f  1 1 


PHOENIX  FIELDSOLVER  WITH  FIXED  JB(R,X) 

PLOT  VERSUS  X 


TIME ,  DATE 
14:04  23 
09/22/81 


I  BEAM  r 
9  993E  +04 

lrise  =■ 

3. OOOE+OI 
ROM  IN  = 
5.000E -01 
ROMAX  = 

1  0001+02 
XLOC  = 

2  OOOE+OI 
XWOIH  - 

I  800E+01 


C*S[  0  >»*•»  13 


-  55  - 


l* 02 


PHOEN  I  X 


F I ELDSOLVER  with  fixed  jb(r,x) 
PLOT  VERSUS  X 


TIME ,  DATE 
M:04:23 
09/22/81 


I  SEAM  — 

9 . 993E+04 
LRISE  = 

3 .000E+01 
ROM  IN  = 
5.000E-0I 
ROMAX  = 
t  000E<-02 
XLOC  i 
2 . OOOE+OI 
XWDTH  = 

t  aooE+oi 


cxst  D  I«»H  u 


X(CM) 


PHOENIX  F I ELDSOLVER  WITH  FIXED  jB(R,X) 

PLOT  VERSUS  X 


T  I  ME  ,  DATE 
14  04  25 
09/22/61 


i  BEAM 
9  993E+04 
LRISE  = 
3.000E+01 
R0M<N  = 

5  OOOE  -01  o 
ROMAX  r  i 

1  OOOE +02  ^ 

XLOC  ^ 

2  OOOE+Ot 

xwotm  ^ 

i  eoot+of 


X(CM) 

case  o  r*m  »5 


-  56  - 


I 


PHOENIX  FIELDSOLVER  WITH  FIXED  JB(R,X) 

PLOT  VERSUS  X 


time,  date 
14:04:23 
09/22/81 


I  BEAM  = 

9 . 993E+04 
LRISC  = 
3.000E+01 
ROM  IN  = 
5.000E-01 
ROMAX  = 

1  - 000E+02 
XLOC  = 

2  000E+01 
XWDfH  = 

1 .80QE+01 


o 

n 

a 

o 


l/T 


X(CM) 


CAS l:  0  f*A>M 


PHOENIX  FIELDSOLVER  WITH  FIXED  JB(R,X) 

PLOT  VERSUS  X 


TIME,  DAU 
14:04:23 
09/22/81 


I  SEAM 

9  993E  *04 
IRISE  = 

3  oooe+oi 

ROM) N  - 
S.OOOf  -01 
ROMAX  r 

r .oooe *02 

XLOC  = 

2  000E+01 
xwoth  - 
1  800E +01 


a 

ii 

a 


X(CM) 

CASE  0  fH/M  17 


57 
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